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We analyse the reconstruction of an unknown pure qubit state. We derive the opti- 
mal guess that can be inferred from any set of measurements on N identical copies 
of the system with the fidelity as a figure of merit. We study in detail the estima- 
tion process with individual von Neumann measurements and demonstrate that 
they are very competitive as compared to (complicated) collective measurements. 
We compute the expressions of the fidelity for large TV and show that individual 
measurement schemes can perform optimally in the asymptotic regime. 



1 Introduction 

State estimation is a fundamental issue in Quantum Information from both 
theoretical and practical points of view. Imagine we are asked to reconstruct 
the unknown state of a quantum system. This can only be done by perform- 
ing measurements on an ensemble of identically prepared systems. With an 
infinite ensemble of copies, the state could be determined exactly. In prac- 
tice, however, we have access to a limited number of copies and the state 
can only be determined approximately^ In this context, three essential ques- 
tions arise: i) what is the optimal measurement scheme?, ii) what is the best 
reconstructed state?, and Hi) how good is the overall estimation process?. 

In recent years a lot of w ork has been devoted to answer these questions 
for different settings! ™ l The optimal strategies, which provide the ulti- 
mate limits that can be achieved, have been identified in several interesting 
cases. However, they usually involve collective measurements (CM), i.e. a 
generalised measurement on all copies at the same time. These, although 
very interesting from the theoretical point of view, are very difficult to imple- 
ment in practice. Far more interesting for experimentalists are individual von 
Neumann measurements for they can be readily performed in a laboratory. In 
this case however, fewer analytical results are knowuPEBI Here, we present 
some theoretical results along these lines. 

We focus on the problem of estimating the most basic quantum state, 
a pure qubit, with physically realizable von Neumann measurements. We 
study quantitatively tomographic inspired schemes )^^ but also consider 
the most general individual measurement procedure, i.e. when depending on 
the previous outcomes, one allows to optimally adapt the measurement on 
the subsequent copyJSI To ease our presentation, we will loosely write local 
measurements for individual von Neumann measurements. Our aim is to in- 
vestigate how good these local measurements are as compared to the collective 



1 



ones. We use the fidelity as the figure of merit quantifying the quality of the 
estimation process and compute the analytical expressions of the average fi- 
delity for large N. Two interesting situations will be analysed which we will 
refer to as 2D and 3D. In 2D the qubit is only known to be a state of the 
equator of the Bloch sphere. This is formally equivalent to phase estimation. 
In 3D no prior knowledge of the qubit is assumed. 

This paper is organised as follows. In next section we obtain the optimal 
state that can be inferred from any set of measurements. This is a general 
result valid for any measurement scheme (either collective or local) and for 
any a priori probability distribution. The 2D and 3D case with fixed (non- 
adaptive) local measurements is studied in section [3J The most general local 
scheme is presented in section We conclude with a summary and outlook 
for further work. 



2 Optimal guess 

The estimation procedure goes as follows. Assume that we are given an en- 
semble of N copies of the qubit state, which we denote by \ri), where n is the 
unique unit vector on the Bloch sphere that satisfies \n)(n\ = (l + n-a)/2 and 
a are the usual Pauli matrices. After performing a set of measurements on 
the N copies of the qubit, one obtains a set of outcomes symbolically denoted 
by x. Based on x, an estimate for \n) is guessed, \M(x)). How well \M(x)) 
approximates the signal state \n) is quantified by the fidelity, defined as the 
overlap 

f n (x)= |(n|M(x))| 2 = ^-U.. (1) 

Eq. Q is a kind of "score": we obtain '1' for a perfect determination (M = ft) 
and '0' for a completely wrong guess [M = —ft). Our aim is to maximize the 
average fidelity, hereafter fidelity in short, over the initial probability and all 
possible outcomes, 

F=(f) = Yl I dnf n (x) Pn (x), (2) 

x 

where p n (%) is the probability of getting outcome x if the signal state was 
\ft), and dn is the a priori probability distribution. For a completely unknown 
qubit, dn is the invariant measure on the two-sphere (on the unit circle in 
2D). Eqs. m and © can be rewritten as F = ±[1 + J2 X V( x ) ' ^{ x % where 



V(x) — J dnnp n (x). (3) 
It is obvious that the choice 

Mix) = (4) 
\V(x)\ 
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maximizes the value of F, which then reads 

F=Ul + Y / \V(x)\^j. (5) 

Eq. (0J gives the best inferred state and Eq. J5J) the maximum fidelity that 
can be obtained for any a priori probability and any measurement scheme 
specified by the conditional probabilites p n (x) . 

In the next sections we show how these simple results can be used to 
improve the estimation procedure. From now on we will only consider the 
2D and 3D isotropic probability distributions given by dn = d(j>/(2ir) and 
dn = sin 6 d6 dtp /(4tt), respectively. A fidelity with no explicit label refers to 
the 3D case. 



3 Fixed local measurements 



Before dealing with local measurements, let us recall some known results of 
the collective schemes. The optimal fidelity for the 2D case iJ^ 
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whereas for the 3D case readsP 

N+l N-^oc _ _j_ 

iV + 2 N 



"CM 



(7) 



These results could, in principle, also be derived from (J5J. Notice that F^^ > 
Fcm VJV, as it should, since in the 2D case we have more a priori information 
about the state than in the 3D case. These results are the absolute upper 
bound for any measurement scheme. 

Let us now turn our attention to local measurements. Any individ- 
ual von Neumann measurement is represented by two projectors 0(±m) = 
(1 ± rh ■ <?)/2, where m is a unit Bloch vector characterizing the measurement 
(in a spin system, e.g., m is the orientation of a Stern-Gerlach). Quantum 
state tomography tells us that, given a large number of copies, von Neumann 
measurements along two (three) fixed orthogonal directions, x, y, (z), are suf- 
ficient to reconstruct the state. 

Consider N = 2"N (3N) copies of the state |n). After INT measurements in 
each axis, we obtain a set of outcomes +1 and —1 with relative frequencies 
on and 1 — on, respectively. This occurs with probability 
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where are the projections of the vector ft in each direction and we have 
used the shorthand notation a = {cti}. Since the expectation value of a is 
(ri\a\n) = n one is driven to propose a guess 

M Ti (a) = , 2a *~ 1 , (9) 

v/E^-i) 2 

where the subscript stands for tomographic. Notice the presence of a normal- 
ization factor such that \Mt\ — 1, therefore Mt always corresponds to a physi- 
cal pure state. Actually, @ is the guess for pure states of maximum likelihood 

procedures The law of large numbers ensures that Mr ""^t^S (n\d\n) = ft, 
but our main goal is to know the rate at which this limit is attained. 

The asymptotic fidelity can essentially be computed by means of the fol- 
lowing systematic approximations (see Bagan et a ;(SlSIl21 f or more details). 
First, use the central limit approximation in (JHJ 

><*NJ ^2^(1 - q) y \ 2q(l-q)J K J 

with q = (1 + rij)/2. Second, transform the discrete sum in into an integral 
using the Euler-McLaurin formula 



1 r 1 



dxf(x) + /(1) - /(0) + //(1) - //(0) -... (11) 



The change of variables r, ; = 2ai — 1, suggested by J5J, proves to be useful to 
simplify the expressions. Finally, use saddle point techniques to evaluate the 
integrals. This just amounts to consider that the value of the integrals is dom- 
inated by the minimum of the exponent in Hl()(l and to expand systematically 
around this point. 

3.1 2D results. 

For the tomographic guess ©, and using the techniques described above, we 
obtain the following asymptotic expression of the fidelity 

4 2D, = l-|i + .-.. (12) 

Note that F^ D ^ approaches unity linearly in 1/N. In this sense, one may 
argue that the tomographic approach is qualitatively similar to the optimal 
collective scheme. Nevertheless, the coefficient of the first correction is a 50% 
larger than the optimal one J§}. 

From our discussion in section[21 we know that there is a better guess that 
can be inferred from the same set of measurements. It is given by (0J, with the 
outcomes labelled by x = {a x ,a y } and the probabilities again given by ©• 
The fidelity is then Fqg = 1/2(1 + J2 a \V{ a )\)- Tne analytical calculation 
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of the large N limit is now more involved, mainly due to the presence of the 
modulus, but can be performed basically with the same techniques.^ It reads 



where OG stands for optimal guess. This is a remarkable result. Provided the 
optimal guess is used, the most basic estimation strategy, namely with local 
and minimal fixed von Neumann measurements, saturates asymptotically the 
optimal CM bound ©. 



3.2 3D results 

The same analysis can be carried out in the 3D case, i.e. when |n) is a com- 
pletely unknown qubit pure state. The calculations are rather more difficult, 
but can be done analytically till the end. For the tomographic guess © we 
obtain 

F T = 1- -— + ..., (14) 
whereas for the optimal guess 

+ (15) 

As expected, Fog > Ft- Notice that, again, the first correction of the fidelity 
goes linearly with 1/N, now with a coefficient very close to one. However, in 
contrast to the 2D case, the improvement of the optimal guess is not sufficient 
to saturate the CM bound for which F = 1 — 1/N + ■ ■ ■. 



4 Optimal local measurements 

The local measurements discussed so far were the most basic ones: fixed and 
minimal. We have not considered yet local schemes in full. In particular, 
classical communication, i.e. the possibility to adapt the orientation of the 
measuring devices depending on previous outcomes, was not exploited. In this 
section, we obtain the optimal scheme in this general setting and show explicit 
results for low N. For large N, we also obtain the asymptotic expression of 
the fidelity. Hereafter only the general case 3D will be considered. 

We need first to introduce a suitable notation to include arbitrary ori- 
entations of the devices and classical communication. Consider the set of 
von Neumann measurements specified by the collection of Bloch vectors 
{fhk}- The set of outcomes x can be expressed as an iV-digit binary number 
x = ijvijv— i • • • *2*i, where ik (= 0, 1). The most general local measurement is 
realized when we allow rnk+i to depend also on the list of previous outcomes 
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ikik—i ' ' "£2*1 = Xk (hence, x — xjv). We thus write m(xk) instead of rhk- 
Note that m(xk) must satisfy the von Neumann condition 

m(lxfe_i) = -m(0xfe_i). (16) 

For any set of outcomes, the optimal guess is given by Q and where now 
the conditional probability is 



N 

Pn(x) = Yl 
k=l 



1 + ft ■ m(xk) 



and the fidelity reads 



F 



/ 2--1 

5 1+ E 



dn n Y\ 



1 + H • m(xk) 



fc=i 



(17) 



(18) 



The optimal scheme is the one that maximizes Ijl8(l over a set of vectors 
{rh(xk)} with the von Neumann constraint (|16J) . 



^.1 Low N cases 

N = 2. Here, there are three independent Bloch vectors vectors: m(0), m(00) 
and m(01) (the other three are obtained using Ea. 116(1 . The first vector m(0) 
is arbitrary and can be fixed at will. The optimal fidelity is then obtained 
by maximizing l(18jl with respect to m(00) and m(01). A straightforward 
calculation yields the following conditions: m(0) • m(00) = = m(0) • m(01). 
Note that m(00) and m(01) do not need to be equal, they are only required to 
be orthogonal to m(0). Substituting back in (|18fl one finds FW = (3 + \/2)/6. 
This is the largest value of fidelity that can be obtained with two copies and 
local measurements. Obviously (3 + v2)/6 < Fcm = 3/4. The optimal guess 
is easily obtained from Eq.Q] M^{x) = [rh{x-i) + m(xi)]/ \/2. This is a very 
gratifying result: M^ 2 \x) is the 'weighted' sum of the outcomes. 

The case N = 3 is very similar. The optimal Bloch vectors, rn(xi), fh(x2), 
rh^xs), are found to be mutually orthogonal. They can be chosen to coincide 
with three fixed (i.e. independent of x) directions. Thus for TV = 3 (as well 
as for N = 2) the optimal estimation schemes based on local measurements 
do not require classical communication. For each outcome x the optimal 
guess is M^\x) = [m{x^) + m(x2) + m(xi)]/-\/3, which is a straightforward 
generalization of M^. The fidelity is F^ = (3 + \/3)/6. These results could 
somehow be anticipated: if 0{rn)\n) ^ we can only be sure that n ^ —m. It 
is then reasonable to explore the plane orthogonal to rh with the next copy of 
\ri). Thus, the optimal Bloch vectors m(xk) tend to be mutually orthogonal. 

The case N — 4 is more complex, since four mutually orthogonal vectors 
cannot fit onto the Bloch sphere . We do not reproduce here the explicit 
expressions of the optimal vectorJ^^. Instead we would like to point out 
some properties of the solution which, in turn, will bring us insight as to how to 
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compute the asymptotic limit. We observe that the optimal Bloch vectors now 
depend on the outcomes of the previous measurements. Therefore classical 
communication does play a crucial role for N > 3. One also sees that the 
third measurement probes the plane orthogonal to the vector one would guess 
from the first two outcomes and analogously does the fourth measurement. 
The fidelity in this case reads = 0.8206, which is just 1.5% lower than the 
absolute CM bound F^ — 5/6 = 0.8333. The maximal fidelities for N — 5,6 
are F^ = 0.8450 and F^ = 0.8637. 



4-2 Asymptotic fidelity 

We can finally compute the asymptotic fidelity of the optimal local scheme. 
Suppose we have performed a (large) number iVo of measurements and ob- 
tained an optimal guess Mq. It is clear that the subsequent guesses will 
hardly differ from Mo- It is also clear from our results of low TV that the fol- 
lowing measurements will basically probe the orthogonal plane of Mq. Hence, 
a good approximation to the optimal local strategy would be to consider: 
a) fixed measurements in the orthogonal plane to Mo (i.e. along two or- 
thonormal vectors u, v of the plane) and 6) a guess of the form M(x) rs 
MqCOSuj + (ucosr + vsmr) smtv, where u> = Xy / (2a u — l) 2 + (2a v — l) 2 , 
tanr = (2a v — l)/(2a u — 1), and A is a tunable parameter. Here a u ,v are 
the relative frequencies of the outcomes as defined in section |3| Note that 
in average to will be small since we expect a u _ v rs 1/2, and only terms up to 
order lo 2 need to be retained. The fidelity can be computed from (0) yielding 

f >i-(i- A) 2 (l-F )-A 2 l ~^~ o Fo) +---, (19) 

where Fq is the optimal fidelity for the first No measurements and the dots 
stand for subleading terms in inverse powers of N and No. If iVo = N@ with 
< /3 < 1,° it is clear that the optimal choice is A = 1, and then F«l - 1/N. 
Therefore local measurements saturate the CM bound at leading order. 



5 Conclusions 



We have obtained the optimal estimation of the a pure qubit state for any 
given set of measurements and and any a priori probability distribution. We 
have focussed on local measurement schemes. For states that are known to 
lay on the equator of the Bloch sphere (2D case), we have explicitly shown 
that, rather surprisingly, the most basic scheme (local and without classical 
communication) saturates the CM bound. This does not happen in the 3D 
case, although the basic scheme yields a fidelity very close to the CM bound. 
We have also obtained the optimal local scheme and shown that indeed the 

a In fact, it can be shown that J3 = 1/2 is the best choice for the partitioning of N. \12\ 
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CM bound is saturated. Furthermore, numerical analysis reveals that the CM 
regime is reached for values of N as low as 12. Our main conclusion is that 
CM do not provide a significant improvement over local measurements. 

Our results can be generalised to other interesting issues, such as estimat- 
ing mixed statespl unknown unitary operations, trace preserving maps, etc. 
For those, the use of local measurements is of outmost interest. 
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